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Abstract 

We study the distribution of transverse energy, Qq, radiated into an arbitrary interjet 
angular region, fi, in high-j»T two-jet events. Using an approximation that emphasizes 
radiation directly from the partons that undergo the hard scattering, we find a distribution 
that can be extrapolated smoothly to = Aqcdi where it vanishes. This method, which 
we apply numerically in a valence quark approximation, provides a class of predictions on 
transverse energy radiated between jets, as a function of jet energy and rapidity, and of 
the choice of the region $7 in which the energy is measured. We discuss the relation of 
our approximation to the radiation from unobserved partons of intermediate energy, whose 
importance was identified by Dasgupta and Salam. 
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1. Introduction 



Jet events at hadron colliders encode QCD time evolution over a large range of length scales, 
from the momentum transfer that produced the jets, to the scale of hadronization. They thus 
afford a window to the transition of quantum chromodynamics from short to long distances, in 
addition to being our most direct probe of new physics at short distances. In this paper, we will 
study transverse energy flow into specified angular regions between high-p^ jets. We will find a 
variety of perturbative predictions, in which soft radiation is sensitive to color and flavor flow at 
short distances. The flow in soft radiation also has an impact on the comparison of measured 
cross sections in a hadronic environment to perturbative calculations. These measurements rely 
on our determination of the energies of jets in the presence of an "underlying event", consisting 
of particles not directly associated with jet production. Because jet cross sections fall so rapidly 
with momentum transfer, they are sensitive to relatively small shifts in energy even at 

the highest pt- This is the case, even though the effects of such shifts are power-suppressed in 
Pt, compared to flnite-order perturbative calculations. An understanding of the underlying event 
requires good control over perturbative bremsstrahlung associated with the hard scattering itself, 
so that the two effects may be separated. 

The influence of the underlying event on jet energy measurements is a complex issue [Q, ^, ^. 
In the absence of further information, it is tempting to think of the particle production underlying 
the jets as similar to that in a minimum bias event, one in which there is no distinguishing hard 
scattering. Energy flow between jets in the flnal state, however, may result from multiple parton 
interactions |]^, |^, which may or may not be the same in jet and minimum bias events, as 
well as from perturbative QCD bremsstrahlung associated with the hard scattering itself [§, 
Marchesini and Webber |TD[ argued that it should be possible to disentangle the underlying event 
from the bremsstrahlung of the hard scattering by studying the flow of transverse energy away 
from the jet direction. They pointed out that the average transverse energy of radiation into 
specifled angular regions is infrared safe. We will develop here a further technique to compute 
the distribution in transverse energy of this radiation. 

The other motivation for this study of interjet radiation is the light it may shed on short- 
distance color and flavor flow, and on the dynamics of hadronization. The relation of the angular 
dependence of interjet particle multiplicity to color flow in the hard scattering has been well- 
established in e"'"e~ annihilation Its use in the analysis of jet events in hadronic scattering, to 



distinguish new physics signals from QCD background, was explored in [jT2|, [T^, and the method 



was verifled experimentally in [Tj], through the analysis of events involving jets and W bosons. 
Here we shall follow [|lOl, and investigate transverse energy flow rather than multiplicity, but the 
spirit of the approach is much the same. 

The discussion below is also closely related to the treatment of rapidity gap ("color singlet 
exchange") events in Refs. |T^, |TB[, for jets at high prp and large rapidity separations. In this 
case, we will study relatively central pairs of jets, with the high-p^ tail of the inclusive jet cross 
section in mind. Correspondingly, we restrict ourselves to a valence approximation. 

We will calculate the distribution for the flow of transverse energy, Qq, radiated into a speciflc 
detector region, fl, away from a set of observed jets and from the beam directions. This is what we 
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mean by interjet radiation.[| The distribution in is perturbatively calculable as a factorized 
cross section, although for ^ Qi with Q the hard momentum transfer, the perturbative 
expansion will in general include contributions of order {1/Qq) a" ln"~ (Qn/Q)- We shall resum 
such terms in a leading logarithmic approximation, in the ratio Pt/Qq, with px the transverse 
momentum of the final-state jets. We shall assume that there are no additional jets of comparable 
energy in the final state. More specifically, we sum inclusively on energy only within cones around 
the beam and outgoing jet directions. We denote by fl the union of all directions outside of these 
cones and outside of Q. To enforce our two-jet criterion, we allow a maximum total energy 
Pt ^ M ^ Qq into Q. More generally, the parameter M can be thought of as shorthand for an 
angular-dependent cutoff on energy flow into Q. 

As shown recently by Dasgupta and Salam [jl7[ |18| in studies of e^e^ annihilation, partons 
emitted into Cl influence energy flow into Q through secondary radiation, producing terms that 
behave as {1/Qq) \n^~^ {M / Qq,) , for n > 2. They call these contributions to the cross section 
"non-global" logarithms. We shall see how they arise in our cross sections, but we shall simplify 
our discussion by neglecting logarithmic corrections in M/ Qq. We comment on the implications of 
this approximation below. The interplay of the definition of jet events and non-global logarithms 
will be the subject of forthcoming work. 

Our analysis is made possible by the incoherence of wide-angle soft gluon radiation from the 
internal evolution of jets, a property of QCD which plays an important role both in proofs of 
factorization |T^, and in the study of jet properties through the coherent branching formalism 



ri| , pm . We will see how non-global logarithms arise in this context from the radiation of gluons 
up to the maximum allowed energy M in region Cl. These hard gluons themselves then act as 
color sources for radiation into Q. Nevertheless, at leading logarithm, because they are emitted 
at wide angles, the eikonal approximation for radiation from an underlying two-to-two hard 
scattering remains. 

We show below how interjet radiation, although independent of the details of the jets, depends 
upon both flavor and color exchanges at the hard scattering, and we will show how to derive the 
distribution of events in terms of Qn for each combination. The choice of Q is free, so long as 
it does not include the observed jets. In this paper, we illustrate our method with the relatively 
simple choice of a region whose area is of order unity in the ?7 — space of rapidity and azimuthal 
angle, outside of the beam-jet scattering plane. 

In the following section, we discuss the class of interjet energy flow cross sections that we 
will study, starting from their collinear-factorized form. We describe the refactorization of hard- 
scattering functions into fully short-distance functions, dependent only on scales at the level of the 
large momentum transfer, ~ p^, and a set of "soft" , but perturbatively calculable, functions, 
which describe evolution at scales intermediate between Q"^ and Qf^. To apply perturbation 
theory, of course, we must assume that Qq is large compared to the scale at which the running 
coupling becomes strong: Qf^ ^ ^qcd- We show how this refactorization leads to renormalization 
group equations for the soft functions, and hence for their Qj^-dependence. These are matrix 



equations describing the mixing of color |2T|. Section 3 deals with the determination of the 



^In fact, our analysis can be applied to a class of kinematic quantities, including the energy itself. We choose 
to work with the transverse energy simply because of its boost invariance. 
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anomalous dimension matrices that control this evolution. We explicitly calculate the matrices 
for initial-state quarks and antiquarks. We show how, at leading logarithm, their calculation 
reduces to the evaluation of integrals over the phase space of radiation into Q. In Section 4, we 
study specific distributions for scales relevant to Tevatron high-p^ two-jet events. We close with 
a brief discussion of our results, and prospects for future progress. 



2. Factorization and Refactorization 
2.1. Dijets and interjet observables 

We will discuss inclusive cross sections involving dijet production in hadronic collisions: 

A{pa) + B{pb) ^ Mpi) + J2{P2) + RniQn) +X^, (1) 

where the jets Jj have high transverse momenta. In principle, the jet axes may be identified by 
any jet-finding convention, but we shall always define the jet energy by flow into some angular 
region, a "cone". In Eq. (|I|), Rn denotes radiation into a region Q outside these jet cones. The 



quantity Qq is a kinematic measure of the radiation into fl. We may, as in Refs. [|T5|, |T6[, take 
Qn to be the energy of particles in Q, but we may also take it to be the transverse energy. As 
above, we refer to the region between Q and these cones by fl. The jets are characterized by their 
transverse momenta and pseudorapidity, rj = ln(cot(6'/2)), where 9 is the angle with respect to 
the beam axis. We will consider two-jet final states for which 

Iml ~ \pt2\ >(5fi > Aqcd • (2) 

We sum inclusively over the radiation Xq in Q, up to a cutoff, M on its total energy in the partonic 
center-of-mass (cm.) frame, with pt M ^ Qq. Thus, the complete interjet radiation into 
fl and fl, although energetic enough to be treated in perturbation theory, is assumed to be soft 
enough that we need not consider the recoil of the high-p^ jets. 



Collinear factorization theorems |jT9| ensure that we may write any such cross section at fixed 
Qq as a sum of convolutions of parton distribution functions (PDFs) that incorporate long- 
distance dynamics, with hard-scattering functions, u, that summarize short- distance dynamics. 



dpi dQi 



^ j dxA dXB <PfA/A (a^A, /^f) 0/b/b {xb, lip) 



d(TAB 

- = 2^ j 

^ IaJb 

XUJf^fg{xAPA,XBPB,Pl,M,QQ,HF,as{lJ,F)) , (3) 



with factorization scale fip. The sum is over parton types, /a,/b- Corrections to this relation 
begin in general with powers of ^qcb/ Qfi^ multiple scattering of partons - part of the 

underlying event referred to above. Our analysis below, which begins with Eq. (|^), is thus 
accurate up to such corrections, and requires the conditions (|^). We take the renormalization 
scale equal to the factorization scale. Because in this paper we will emphasize ratios of cross 
sections rather than their size, fixing /iren = /^f = Pt is meant to be taken as an indication of the 
order of these scales, rather than the promotion of a specific choice. 
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2.2. Refactorization of the partonic cross section 

As stated above, we define the energies of liigli-pr jets by tlie energy flowing into specified cones, 
while the region Q in which we measure Qq is outside these cones. Soft gluon emission outside 
the angular extent of the jets decouples from the kinematics of the hard scattering, and from 
the internal evolution of the jets [0, As a result, we may express ujf^f^ in Eq. (|^) as a 



sum of terms, each characterized by a definite number of jets produced at the hard scattering. 
To lowest order in a^, only two jets are possible. Thus, at leading order, the production of the 
high-pT jets is given by the set of Born-level 2^2 partonic processes, which we label f, 

f : fA + fB^fl + f2. (4) 

We distinguish qq qq (/i = q) from qq — > qq (/i = q). We may now write the single-jet 
inclusive cross section at fixed as 



with fj = T] — {l/2)\n{xA/xB) the jet rapidity in the partonic center-of-mass, and s = xaXbs. 
We neglect the effects of recoil of the observed jet against relatively soft radiation. The short- 
distance function da^^^ /dfjdQfi is an expansion in asipr) that begins with the lowest order (LO) 
cross section, 

dT]dQn dr] 

but which includes potentially large corrections associated with logarithms of Qn/pr- In the 
two-jet approximation described above, we work at leading order in the hard scattering, with the 
factorization scale fip = Pt- The sum runs over all possible partonic subprocesses f, defined as 
in (|). We emphasize that corrections to Eq. begin at 0(as(pr)), due to processes with more 
than two energetic jets in the final state. Taking these into account would require generalizing 
the set of flavor flows f in Eq. to 2 — >^ 3 and beyond. Such a development will be informative, 
but we shall not pursue it here. 

The hard-scattering function for the subprocess f may now be refactorized [^, ^ , to exhibit 



the interrelation between the partonic hard scattering and soft radiation outside all jet cones. 
For this purpose, it is convenient to treat the partonic and hadronic cross sections integrated 
over the transverse energy radiated into region fl, up to a fixed value Qn, 

rfa(^) (Pt, V, f^F, M, Qn, asifip)) _ f^^ d&^^^ {pr, V, fJ^F, M, Q'^, ^^(/iF)) 



di] Jo dfjdQ'n 

dcJAB {pt, V, f^F, M, Qn, as{fiF)) 



drjdpj 



J dXAdXB (f)fA/A{xA, IJ'f) <PfB/B{xB, fJ'F) (7) 



\ rfo-ffl {pt, r), fiF, M, Qn, asjl^F)) 
2 cosh fj I dfj 



5 



As we shall show, the relevant refactorization of the hard scattering function of Eq. (|^ is 



dr] ^ \ fi fi J 

(f) 

where Hj£ is a matrix of short- distance functions, independent of Qq, while all of the Qq- 
dependence is factored into dimensionless soft functions S\^j. Here yU is a new refactorization 
scale, chosen between px and Qq, with dynamics on scales below fi incorporated into S^^\ The 
soft function is independent of fip, which separates da/dfj from the parton distribution functions. 
The functions S}^j, which summarize dynamics of partons with energies at the scale Qq, depend in 
general on the scalar products Pi-Pj between the lightlike 4- velocities of the incoming and outgoing 
energetic partons in the hard scattering, which radiate soft gluons. If we fix the scattering plane 
such that jet 1 is located at = 0, the kinematics of two-to-two scattering reduces this dependence 
to a single variable, which we may take as fj, the rapidity of parton 1 (see Eq. (§)) in the partonic 
cm. frame. Note that S\^j is also a function of the maximum energy, M, radiated into the 
complementary interjet region, Cl. 

In Eq. (H), the hard-scattering functions Hjl begin at order a^, while the soft functions 
Sli begin at zeroth order. To leading logarithm, we need only these lowest-order contributions, 
denoted below as H^^'^^ and S^^'^\ respectively. The combination of these approximations is 
related to the lowest-order cross section by Eq. (||), 

= E Hf, iPT, V, asifip)) S^i'^ . (9) 

L,I 

Recalling that we expect corrections to Eq. to begin at NLO in asipr) for inclusive jet pro- 
duction, Eq. (ID should hold to a corresponding accuracy in the hard-scattering function Hjl- 
Additional corrections to the refactorization (||) may be due to multiple scattering. These correc- 
tions are suppressed by powers of ^qcd /Qh- large Qn, therefore, bremsstrahlung dominates 
the inclusive cross section. At low Qq, the NLO differential cross section for bremsstrahlung 
diverges as l/Qn- As we shall see, however, the resummed integrated cross section, Eq. (|^), is 
better behaved, and may be extrapolated smoothly as far down as Qq = Aqcd, where it van- 
ishes. In these terms, the perturbative component of interjet radiation, which dominates for hard 
radiation, joins smoothly to the region of low radiation, and may thus be meaningfully compared 
to data [ITU . 



2.3. The soft function 

The proof of Eq. (^ at the desired accuracy, leading logarithm in Qq, follows standard arguments 
in factorization [p!9| , p3| . An essential ingredient in proofs of collinear factorization is that wide- 
angle soft radiation decouples from jet evolution [^. Soft radiation away from jet directions is 
equally well described by radiation from a set of path-ordered exponentials - nonabelian phase 
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operators or Wilson lines - that replace each of the partons involved in the hard scattering (four, 
in our case), 



$J/)(oo,0;x) 



I { 

Pexpi-ig / d\ P-A^f^ {X(3 + x] 



Pexp(-ig r d\(3'-A^f'\\l3' + x] 



(10) 



where P denotes path ordering. The first line describes an outgoing, and the second an incoming, 
parton, whose flavors and four-velocities are / and /5, and /' and respectively. The vector 
potentials A'^^'^ are in the color representation appropriate to flavor /, and similarly for A^^ \ 

Because wide-angle, soft radiation is independent of the internal jet evolution, products of 
nonabelian phase operators, linked at the hard scattering by a tensor in the space of color indices 
generate the same wide-angle radiation as the full jets. The general form for these operators, 
exhibiting their color indices, is 



{d^} 



X c 



.(f) 



<l>!i-H0,-oo;x),„,, $)/;^(0,-oo;x),,,,, . 



d2,di;dA,dB 



ifB), 



The c/ are color tensors in a convenient basis. Examples will be given below. The perturbative 
expansions for these operators are in terms of standard eikonal vertices and propagators, and 
have been given in detail in Refs. pT| , |22| . In these terms, we define 5*^^^ as [| 



'-'LI 



r],n, — , — ,as(/i) 



dQ'n E E ^(Q 



in) 

n 



x(0|f 



t 



n)(n|T[iyf (0)|J |0). (12) 



The sum is over all states n whose transverse energy flow into Q is restricted to equal Qq, and 
the energy into fl is less than or equal to M. The indices L and / refer to the color exchange 
at the hard scattering between the partons in reaction f, as built into the definitions of the W^s, 
Eq. (|TT]). The matrix elements in Eq. (|r^) require renormalization, and we may identify the 
corresponding renormalization scale with the refactorization scale of Eq. (|]). 

The Qn-dependence of the matrix S^j in Eq. (^2|) is the same as in the full partonic cross 
section, up to corrections due to differences between the jets and the nonabelian phase operators. 



Radiation within the jets, however, is treated inclusively in the sum over final states in (|T2D. As 
observed above, radiation outside the jet cone is emitted coherently by the entire jet JlT], |l^, ^ 



Any logarithmic enhancements associated with (angular ordered) radiation within the jets, or 
with radiation collinear to the eikonal lines, are thus independent of Qn, and cancel in the 
sum over states. This cancellation occurs at fixed momenta for all on-shell particles within the 

^This definition of S^^^ differs slightly from that of Rcf. H, ||, in that it is an integral over Qo, rather 
than differential in Qq. 



7 



jets p5[. As a result, in the leading logarithmic approximation, differences between the 



eikonal and exact cross sections are given by Q^-independent expansions in which may 

be absorbed into the hard functions Hlj. The iJ's themselves are thus the coefficient functions 
that match the effective eikonal operators Wj to the full theory. Therefore, Eq. @, with the 
soft function defined as in Eq. (|T2|), reproduces the Qo-dependence of the original cross section, 
to leading logarithm, the level of the original collinear factorization form restricted to a 2 ^ 2 
hard scattering, Eq. (|^). 

We shall see in the next subsection that the color flows associated with 2^3 and beyond 
play an important role in the analysis of the soft matrix S^^j itself, even at leading logarithm |1T7| . 
These effects do not, however, affect the basic refactorization of eikonal from partonic degrees 
of freedom, Eq. (H). In addition, as noted above, beyond leading logarithm the sums over / and 
L in the factorized hard scattering cross section, Eq. (^, must also be generalized to include 
the possibility of 2 ^ 3 and other scattering processes, in which more lines emerge from the 
short-distance processes. To include these higher-order processes, we would generalize the 2^2 

/r\ I— K 

eikonal cross sections of S\^j of Eq. (|12]) to 2 — > 3, by including three (or more) outgoing eikonal 
lines in W^^\ Eq. ([TTp. 



2.4. Resummation for the soft function 



The partonic cross section a must be independent of the refactorization scale, fi in Eq. (§), at 
which the soft (but still perturbative) wide-angle radiation, sensitive only to overall color flow, 
and the truly short-distance dynamics are separated: 



d 



dfj 



0. 



(13) 



This condition applied to the right-hand side of (|^) leads to a renormalization group equation, 



S 



JI 



(14) 



with, as usual, Qs = i/ivraj. Here (^r5''(r))^ may be thought of as an anomalous dimension 
matrix |2^. These anomalous dimension matrices, and their eigenvalues, depend on the cm. 
rapidity, fj, of the partons that participate in the hard scattering. Equivalently, they depend 
upon 6, the scattering angle in the partonic center-of-mass. We always fix the observed jet at 
= 0. 

Because the soft function depends on the cutoff M for radiation into 0, the Qn-dependence 
is related to the anomalous dimensions by 



d_ 



JI 



M 



95« 



LI 



DM 



(15) 



To make this equation explicit, we must interpret the derivative with respect to the cutoff. 
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The derivative with respect to M fixes the energy fiow into f2 equal to its maximum value, as 



may be seen from Eq. (0). This ensures that at least one gluon with energy much larger than 
appears in VL. Because the soft function begins at order the result may always be written 

as 

= - (A4"(^.f!))*^sS - SS (A4"(f,,fi))^^ +5 (f,,^l.^.a,(^)) . (16) 

in terms of an effective anomalous dimension matrix AF and a remainder term, S. Both AF and 
S are computed from real-gluon emission, and hence depend on the geometry of the region VL in 
which soft radiation is measured. The remainder, S begins at order a^, and in general includes 
contributions proportional to In • The logarithmic enhancements in these contributions 

are associated with the ratio M/Qq. In these contributions, energetic gluons, of energy ~ M 
produce secondary radiation that appears in VL. In iS, this radiation resolves the color flow of the 
energetic radiation from that of the original eikonal lines that deflne the soft function in Eq. (p^). 
As such, we cannot absorb these contributions into an anomalous dimension matrix for 2^2 
scattering, but must consider more general color flows. To preserve the simplicity of a single 
matrix anomalous dimension, we shall simply neglect S below. In a sense, this is equivalent to 
neglecting logarithms of M/Qq^ compared to those in px/Qn- This is by no means a systematic 
approximation, however, and we must regard the resulting formalism as a somewhat simplifled 



model for the actual cross section. The study of Ref. |T8| for e"'"e~ annihilation, however, suggests 
that the effect of neglecting S does not change the behavior of the cross section qualitatively. 
These reservations notwithstanding, we now deflne 

(r?(^,^^)),, = (r?(^)),, + [^^f{v.^)),, , (17) 

which depends on the geometry of the region VL in which soft radiation is observed. Then, 
neglecting iS, we have an equation for the Qn dependence of the soft matrix 

(q.^^ - I3{a.)-Q s'l] = (r<"(f,, a))[^ sf, + sS (r?(^. a))^^ . (is) 

As discussed above, Eq. (|^) holds at leading logarithm in /i, and hence Qn/l^', consistent with 
the accuracy of the underlying factorizations, Eqs. @ and (H), but neglects corrections at the 
same power in ln(M/Qn). In the following, we therefore drop the argument M in the cross 
section and soft function. 

To solve Eq. (|l^), we go to a basis for the color matrices c^P that diagonalizes TP{'q,Vt), 
through a transformation matrix i?. 



where 



Xf{f,,^l) = '^\f'\f,,^l) + ... (20) 



TT 
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are the eigenvalues of r^^''(j7, Here and below, Greek indices /3, 7 indicate that a matrix is 
evaluated in the basis where the soft anomalous dimension has been diagonalized. Thus, for the 
soft and short-distance functions we also write. 



oil) 



(f) 



7/3 



S 



(f) 
LK 



L/3 



-1 



(21) 



The transformation matrix R^^^ ^ is given by the eigenvectors of the anomalous dimension matrix. 



K 



(22) 



The solution to Eq. (p^), which resums leading logarithms of Qn/ f^, is straightforward. We 



introduce a combination of eigenvalues of Tg , E^^l^fj, Q), given by 



(23) 



in terms of the lowest-order coefficient of the beta function, /3o = (llA^c — 2nj)/3, for Nc colors 
and Uf quark flavors. The soft function is then 



S. 



(f) 

7/9 



17, fi, — ,a^(/i) 



{v,^A,as{,Qn)) exp 



5'5](r),fi,l,a,(Qf,)) 



J U 



27r 



(24) 



where in the second form we have reexpressed the result in terms of one-loop running couplings. 
We now set the refactorization scale, /x, equal to the transverse momentum of the observed jet, 
Pt, and find, for the partonic cross section. 



da^^^ {Pt, V, fiF, Qa, «^(/^f)) 
dfj 



A 7 

xS^'l{fj,n,l,as{Qn)) 



'Osipj 



(25) 



At leading logarithm, we approximate if'-^^ by if^^'^^ and S*^^-* by S^^'^\ as observed above. Then, 
for = this expression reduces to the Born cross section for inclusive jet production, as in 
Eq. (H). The corresponding differential cross section in Qq, valid to leading logarithm, is 



Qn 



da^^^ {Pt, 17, yUF, Qn, asi^r)) 



dfjdQi 



/3,7 



2n 



as{p 



'T 



"s(<5n) 



(26) 
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For consistency, we must recover at lowest order the eikonal approximation to the complete 2 — 3 
cross section at fixed from (^61) . We shall see how this result emerges in the following section, 



where we describe the calculation of the exponents 

Once the matrix R^^^~^, Eq. (p2D, is determined in a specific basis for color exchange, the soft 
and short-distance matrices in Eq. ( p6D are found by applying Eq. (pT]). The lowest-order soft 
matrices S^^^ reduce to color traces only, while at leading logarithm, the hard matrices H^xl^ 
are given by squares of the corresponding leading-order partonic 2 — > 2 amplitudes, projected 
onto the original color basis. 

We are now ready to introduce a convenient measure of soft radiation into region Q in our 
dijet events. Consider the ratio 

= — , (27) 

where the numerator is the inclusive hadronic cross section given in Eq. (0) above, and the 
denominator is the lowest-order dijet cross section for the same hadrons. 



drjdpT 



J dXAdXB (pfA/A {^A, /Uf) (pfB/B {xb, I^f) 



7,L 



Here we have used Eq. (^) to relate the LO hard-scattering function a, in collinear factorization, 
to the lowest-order, refactorized hard and soft functions, H and 5*. By comparing this expression 
with Eqs. (0) and (^) for the numerator in (|27|), we see that, so long as Re E^ij > 0, pab 
vanishes for small Qq (strictly speaking, at = Aqcd), increases monotonically with Qq, 
and reaches unity at Qn = pt- This behavior, due to resummation, is to be contrasted with 
the corresponding NLO expression for the same quantity, which diverges toward minus infinity 
for vanishing Qq, where the virtual contribution dominates. The quantity pab predicts the 
distribution of bremsstrahlung into region Q. 



3. The Anomalous Dimension Matrices 

We now turn to the calculation of the anomalous dimensions Tg\ introduced in Eq. (|1^) above. 
Equation (|1^), in turn, followed directly from the refactorization of the hard-scattering function, 
Eq. (^, with an extra contribution from real-gluon emission into Cl, as in Eqs. (|1^) and (0). The 
refactorization scale, p of Eq. (|^) was identified with the renormalization scale for the operators 
W in the definition of the soft functions S^^\ Eq. (|I2p. In this section, we shall derive these 
anomalous dimensions from the renormalization of the soft function. To compute the anomalous 
dimensions at one loop, we may identify the cutoff M for real gluon emission in the definition of 
5*^^^ with the refactorization scale, p. We thus suppress the additional argument M in for this 
calculation. We will work in Feynman gauge. 
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The renormalization of multi-eikonal vertices has been discussed in some detail in Ref. ET 



and we shall follow the method outlined there, with an important difference. In the anoma- 



lous dimensions Tg were computed in axial gauge, after dividing an eikonal scattering amplitude 
by eikonal self-energy functions. This extra factorization eliminated double poles in dimensional 
regularization, which are associated with coUinear emission by the nonabelian phase operators. 
In axial gauge, these singularities appear only in self-energy diagrams. In the present case, such 
an extra factorization is not necessary, because S'}f], defined as in ([T^) is free of coUinear singu- 
larities associated with the eikonal lines, after the sum over intermediate states. The remaining 



ultraviolet divergences in S^^j may be compensated through local counterterms, just as in |£1 



From these counterterms, we may simply read off the entries of the anomalous dimension matrix 
in an MS renormalization scheme. 

Given the definition of the soft function, Eq. (|l5), the counterterms for S in Eq. (H) are 
calculated only after a sum over final states. They therefore depend on the choice of fl, and the 
scattering angle of the underlying process. The following two subsections describe the one-loop 
calculation of S\^j in this approach. 



3.1. Renormalization and color mixing 

To begin, we consider the diagrams shown in Fig. Figs. |I|a and |T]c are virtual corrections, 
while |p3 shows the corresponding diagrams for real gluon emission. To compute the counterterms 
associated with these diagrams, we fix the transverse energy flow = in Eq. 

As is characteristic of eikonal diagrams. Figs. p]a-c are all proportional to scaleless integrals, 
which vanish in dimensional regularization. At Qn = 0, therefore, S);^} is given entirely by its 
one-loop counterterms, which we denote as 

-Slf^ (4'i)(r),fi,a.(/i),e))^^, (29) 

where S^^^ is the zeroth order matrix, which, by Eq. ( |12|) is a set of color factors, independent 
of 17 and fl. The elements of Z are computed in D = 4 — e dimensions. It is important to note 
that the e poles, given by the one-loop counterterms, have the interpretation of infrared, rather 
than ultraviolet, divergences in S^^j, and that they cancel infrared divergences from real-gluon 
emission into region Q.^ Following the standard analysis, the one-loop anomalous dimensions are 
given by 

(Tfir),^)),j = -«.^ Res,^o (4'''^ (r/, fi, a.(/i), e))^_^ . (30) 

The calculation of the Z's is thus essentially equivalent to the calculation of the anomalous 
dimensions. To carry out these calculations, however, we must specify a basis of color tensors, 

Cj. 

■^More specifically, the divergences from virtual diagrams are proportional to 5{Q'^) in ([l^), and cancel real- 
gluon emission as integrable distributions. 
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Figure 1: Diagrams for the calculation of the anomalous dimension matrix, as in Eq. (^31). The 
double lines represent eikonal propagators, linked by vertices c) and , in the amplitude and 
its complex conjugate. The vertical line represents the final state. 



For the two-to-two partonic processes, Eq. (^, we denote by the color index associated 
with parton (or eikonal line) i = A,B, 1, 2. As above, capital letters, I, J . . . are indices in the 
space of color exchanges, spanned by a chosen basis for each partonic process. We will restrict 
the discussion here to the calculation of anomalous dimensions when the incoming partons are 
quarks or antiquarks. Convenient bases for describing the color flow for quark processes are 
i-channel singlet-octet bases, given by 



for qq — > qq, and 



Cl 
C2 



2N, 



(31) 



Cl 
C2 



for qq qq. The process qq 



^rA,ri^rB,r2i 



1 



rA,T2^rB,ri • 



2j^^-'A,..-.B,.. ■ 2^ 
gg is best described in terms of the color basis 



(32) 



Cl 
C2 
C3 



^ri,r2 ^rA,rB ' 
dr2ric (^F)-rR» 
'^fr2r\c (Tp 



TBTA 



(33) 



'TBTA 
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where Ci is the s-channel singlet tensor, and and C3 are the symmetric and antisymmetric octet 
tensors, respectively. Bases for other possible partonic processes with gluon eikonal lines have 
been given, for example, in Refs. |T6| , |21| . 

Let Z^"^^^ denote the contribution to the counterterms from all the one-loop graphs in which 
the gluon connects eikonal lines % and j. In this notation, the calculation of Fig. |I| gives us 
We emphasize that the Z's, and hence the anomalous dimensions, are to be calculated only 
after combining real-gluon emission diagrams with virtual diagrams. The Z's are constructed to 
give local counterterms that cancel only those ultraviolet divergences that are left over after the 
real- virtual cancellation has been carried out. 

To find the z\^p^s in the color bases (^) and (|32|), we rewrite the various one-loop virtual 
diagrams in terms of the original color basis, using the identity shown in Fig. ^, 



W^i = I {SaS.k - ^S.Ai ) (34) 



for quark processes. For scattering processes involving gluons, many useful identities can be 
found, for example, in p6[ . This procedure results in a 2 x 2 matrix decomposition for scattering 



involving only quark and antiquark eikonal lines, describing the mixing under renormalization of 
their color exchanges. The annihilation of a pair of quark and antiquark eikonals into two gluon 
eikonal lines gives a 3 x 3 matrix structure, while for incoming and outgoing gluonic eikonals, we 
get an 8 X 8 matrix pT| . 

For a given Z^^^\ the momentum-space integral appears as an overall factor. It is then 
convenient to introduce the notation 

{Zl'^) = (35) 

where the (li are the coefficients that result from the color decomposition of the virtual diagram, 
and where the to^'^^^s include the ultraviolet pole part of the momentum space integral for the 
soft function 5* and remaining overall constants. Defined in this fashion, the sign of each u^^^^ 
depends on the fiow of fiavor in the underlying process (see below). From ( pOD and (^) the 
relation between the cu's and the anomalous dimension matrices is 



- E cir ^^'^^ ■ (36) 

{(ij)} 



As an example, we give the color decomposition of F*^"^^) in qq qq, for which we find, using 
Eq. (H), 

C(^^) = ( ° i|) . (37) 

The step to the anomalous dimensions is trivial, following Eqs. ( pOD and (pGl). We define 

r^'^'> = -euj^'^\ (38) 
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= 1/2 




1/(2N,) 



Figure 2: Color identity corresponding to Eq. 



Then we may summarize the full anomalous dimension matrices ( |36D in terms of the following 
combinations, 

^(f) ^ p{Al,f) _^ p(B2,f)^ 



7 



(f) 



p(A2,f) _^ p(Bl,f)_ 



In the basis (^Tj), the anomalous dimension matrix for the process qq — >■ qq is given by 



^(qq^qq) _|_ ^(qq^qq) 



(39) 



(40) 



Here and below, we have suppressed the arguments of the momentum factors, and 7 for 
compactness, but they all depend on 17 and on Vt. The color decomposition for the process qq —>■ qq 
is the same as above, Eq. (140|). The momentum parts change, however: p(i^^^i') = —^im^w) ^nd 
ryiqq^qq) = _p{qq^qq)_ 

Similarly, for qq 



r 



(qq-^qq) 
s 



qq in the basis (|3^), we have 

(J^piqq^qq) ^ (a^ll^'^l) + ^(ll^^q)^ 



(yiqq^qq) _|_ ^{qq^qq) (j^^iqq^qq) 3_ + _|_ Aqq^qq) 



(41) 



|(/? + 7) \ 
f (/5 + 7) 



(42) 



Note that certain momentum factors for quark-quark and quark-antiquark differ by signs, due 
to the eikonal Feynman rules (see below. Sec. 3.2). 
For qq gg we find in the basis ( ^3)) , 

where ^ = ^ (/5 — 7) — 2^r*^^^) + ^F^^^^ Here we have suppressed the superscripts, and define 
the F*^*-'^ as well as a, (3 and 7 for qq — > gg to be the same as those above for qq — > qq. We 
have done so because, for processes involving gluons, the factorization into momentum and color 
factors as in (|35|) is ambiguous. The gluon-eikonal gluon vertex has both momentum and color 
parts, and the sign of each depends on the orientation with which the exchanged gluon attaches 
to the gluon eikonal line. Their product, and the anomalous dimension matrix, of course, is 
independent of the way the graphs are drawn. 
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3.2. Momentum integrals 

We now compute the momentum integrals cu^^^^ that contribute to the matrix of renormahzation 
constants (^Z^^^ at one loop, through Eq. (|35D . The momentum factor for the virtual diagram 
in which the gluon line connects eikonal lines i and j, such as Fig. |l|a, may be written as |22 

4«' m. A.. A,) ^ / ^ ^ + • («) 

where here and below the subscript P.P. indicates that the integral is defined by its ultraviolet 



pole part. In (43), 6i = 1(— 1) for momentum k flowing in the same (opposite) direction as the 
momentum flow of line i, and Aj = 1(— 1) for i a quark (antiquark) line, while, as noted above, 
we define the momentum integrals for gluons to have the same signs as for quark-antiquark 
scattering. For example. Fig. |I|a for qq —>■ qq has i = A, j = 2, with = 1, A2 = —1 and with 
= ^2- In this notation, each of the real gluon emission diagrams, such as in Fig. |l]b, takes the 
form 

— * 

with Si and Sj defined by the corresponding virtual diagram. The function Q(k) = 1 when the 
vector k is directed into region Q, and is zero otherwise. Both real-emission diagrams give the 
same answer, while the two virtual diagrams are complex conjugates of each other. 

The sum of the virtual diagram in the amplitude with either of the real diagrams gives: 



d 



-9s 



p.p. 

+AiA,6i6j ^ ^ (1 - 



(45) 

t 

In the last line, we have recalled the definition (|38D. The term with the fc-integral is the real 
remainder of the cancellation of the real and virtual diagrams, while the imaginary part is nonzero 
only for (ij) = (12) or (AB). At Qq = 0, the real diagrams oj^i^ get contributions only from 
gluon emission outside of fi, and cancel contributions from Uy''^ as shown in Eq. (^5]). The 
counterterms are identified after this cancellation has been carried out, and are implemented 
as local subtractions at the color vertices, cj and c^. This requires that the counterterms, and 
therefore the anomalous dimensions, be complex. 

Given the anomalous dimensions computed in this fashion, and the lowest-order hard and soft 
functions, we find the differential cross section for interjet radiation into region Q from Eq. (p^). 
Note that for Qq 7^ 0, the cross section may be fixed uniquely at one loop by the condition that 
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the integral of the real part down to Qq = must cancel the divergences of the virtual diagrams, 
that is, the counterterms identified above. We recognize that the result is consistent with Eq. 
([26| ) expanded to one loop, after a transformation to the diagonal basis for color exchange. The 
inclusive cross sections are given similarly, by Eq. (|25|) . 



4. Out-of-plane Radiation 

As an illustration of our method, we introduce a region Q that is fixed out of the scattering 
plane, as a "rectangular patch" in the space of rapidity y and azimuthal angle 0. Other choices 
are possible of course, and will be explored elsewhere. 

Taking the observed jet at = 0, and the recoil jet at = vr, we define Q in the hadronic 
center-of-mass frame by 

(pmin ^ ^ (j^max i 
Vmin < y <r] 

max ; 

(46) 

where the boundaries are chosen such that the patch is located outside the incoming beams 
and outgoing jets. The distribution of energy radiated into the out-of-plane patch is determined 
by the eigenvectors and eigenvalues of the anomalous dimension matrices, Tg\ Eqs. (|40|) , (W^), 



and p^ ) above. Given the definition ( |4^ ) for Q, and the expressions ( |45D for the diagrammatic 
contributions, we simply need to evaluate the one loop integrals over the loop momenta restricted 
to this phase space. 



4.1. Entries for the anomalous dimensions 

For the patch region of Eq. (|46|) , all the phase space integrals of Eq. (^), and hence all the entries 
of the anomalous dimensions, F^. , in Eqs. (^)-(^) can be expressed in terms of the following 
combinations of logarithmic and dilogarithmic functions. 



I{x, y, z) = arctan(z) In 




We note that this function has a finite limit to z = ±i, where singularities of the dilogarithms 
cancel those of the arctangent: 



I{x,y,±i) = 




For the anomalous dimension matrices relevant in the valence quark approximation, we need 
explicit expressions for F(^^'^) and F^^^.f)^ addition to a^, p^^^ and 7^^), defined in Eq. 
With variables in (^6D transformed to the partonic center-of-mass frame, the relevant arguments 
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of the J's in Eq. are: 



Xmin = tanh + Vrmn ^ tan{(j)max / 2) , 

tanh tan(0ma^/2), 
coth — tan(0maa;/2), 

I/max = coth tan{(f).max/'^), 

z = tan 0^„,/2 . 49 

cosh (17) 



In these terms, we have 

271 



'1 

^2 



^ \'lmax I Imin ) yr max YminJ 

as 
27r2 



as 



p{ Ql Ql) ^ I (^Xmim •^niax J ^) I {p^mim •^max ^ Z) I i^Xminj •^maxy Z ) } {(pmax ^ (pmin} 

0) + I{y 

min 5 Vmax 5 

z) + I{y 

mini y max 1 ■^^ ) } {4^ max 4^ min} 

- 2nH, 

27]-2 

(3^^^ ^ {^( ymini ymaxi 0) 2 /( ymini ymaxi t'^-'^(*^max/2) ) } {4max ^ 4min} 

Ois 

+ {I {^XfYiini -^maxy 0) 2 I{Xmini ■^maxi tan(0maa; /2) ) ]■ {(pmax 4min\ ; (^^) 

while 7(99^'?^) is found from jS^i'i^i'i) through the transformation, 

ry{l~ '^'^^ (^fj,^ifi J fjmaxi (pmini (pmaxi ^i) ( Vmini Vmaxj (pmini 4maxi ''}) • (^■^) 



Expressions for quark-quark scattering differ only by overall signs, determined by the factors Aj 
in Eq. (||). 

The anomalous dimension matrices are found by substituting these results into Eq. (0) for qq 
scattering into qq, or (^) for their annihilation into gluons, or ( ^T]) for the quark-quark case. It 
is then straightforward to find the corresponding eigenvalues, \^^\ Eq. (20), and transformation 
matrices R^^^-\ Eq. (|2|). The resuhs are somewhat cumbersome, and we do not present them 
here. In Figs. § and we plot the real and imaginary parts of the exponents Eap, Eq. (^), 
for qq qq and qq — » gg, as functions of f), the rapidity of the observed jet in the partonic 
center-of-mass, for a typical out-of-plane patch Q, Eq. (|46|). For these plots we have chosen a 
configuration with xa = a^s in the factorization, Eq. (|^), so that the partonic and hadronic cm. 
frames coincide, with rjmax = —Vmin = 1, (l>min = 7r/4 and ipmax = 37r/4. Note that Eap changes 
with xa and as f2, fixed in the hadronic cm., is boosted in the partonic cm. The real parts of 



18 





.0 1.5 



Figure 3: Real and imaginary parts of Eaf^ for qq qq with Q chosen as in Eq. (^), with 
4>min = 7r/4, (j)max = 37r/4, and rimax = -Vmin = 1- For the real parts we plot only the diagonal 
elements, Re Eaa, since the interference terms, Re Ea/3 with a ^ P, are just the averages of 
the former. Note Im E^a = 0, and that for an n x n anomalous dimension matrix, only n — 1 
imaginary parts of the E^p with a ^ (3 are independent. 



Ectp are all of order unity in the range of i) shown. We expect them to be stable against higher- 
order corrections, so long as the size of the patch is of order unity. For very small patches, we 
would in general anticipate large logarithms associated with an incomplete cancellation between 
real and virtual corrections. 

The only remaining ingredients necessary to compute the resummed cross sections, Eqs. ( pSj) 
and (|26|) are the hard and soft functions at lowest order, to which we now turn. 



4.2. Hard and soft matrices at LO 



Explicit hard and soft matrices in color space at LO for the partonic processes qq — » qq and 
qq —* qq have been given in Ref. |TB[. We exhibit them here, along with those for qq gg, with 
a trivial change in overall normalization relative to [0. To use them in the resummed cross 
sections, Eqs. (|^) or (^), it is necessary to apply the transformation rules (^), to take the 
results quoted below to the ^-dependent diagonal basis. 



4.2.1. Hard and soft matrices for qq qq at LO 



In the basis ( pTD the decomposition of the Born level hard scattering for same-flavor qq 
color space is given by the matrix 



( (QeS"^^ Of. 



qq m 



(52) 
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Figure 4: Real and imaginary parts of Eap for qq gg. 



where Xi, X2, and Xa are defined by 



Xi 



52 



X2 = A/^cTT 

_ S^ + U^ 1 t2^M2 2 m2 

- ^^^Nl^^'KJi- ^^^^ 

The step to the unequal-flavor cases, qq' — >• qq' and — >• g'g' is straightforward, by dropping the 
s-channel terms for the former, and the t-channel contributions for the latter. The matrix for 
qq — >■ qq is found from (|52D by exchanging t and u. The lowest-order soft matrix is just given by 
the trace of the color basis, Sf} = Tr {c^cj), resulting in 

4.2.2. Hard and soft matrices for qq qq at LO 

iFoT this process we get, in the basis Eq. (^), a hard matrix related to the one for qq — > qq, Eq. 
([5^), by the crossing transformation s ^ u. The functions Xi? X2 and Xs here are given by 



Xi 



s^ + i^ 

V? 

S2 s2 + P 



X2 = N,^- 



tu 

_ s^ + u" 1 s^ + i^ 2 52 

^2 IS^ MctU 
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Again, qq' — > qq' is found by keeping only the t-channel terms. The lowest-order soft matrix is 
the same as for qq qq. 

4.2.3. Hard and soft matrices for qq gg at LO 

In the basis (^3|), the Born level hard-scattering matrix is 

/ mXi ifcxi 1^X2 \ 



Xi 



\ JV7^2 X2 



where xi, X2, and xs ai's now given by 

Xi = 

X2 = 
X3 = 

The zeroth order soft matrix is found to be 



iu 



2Nl 
Nl-A 








X2 
X3 ) 



(56) 



(57) 



(5^ 



4.3. Cross sections 

Given the hard and soft functions, and the entries of the anomalous dimension matrices, we 
have all the ingredients necessary to calculate the inclusive cross section (P3[), its differential 
counterpart (^), and the ratio p, Eq. ([27|), for pp scattering in valence approximation. For a 
given choice of f2, such as the patch described above in terms of azimuthal angles and rapidities, 
we find first the functions r*^*-''^^ (5^^^ and 7^^^ for (f) = (gg qq), {qq qq), {qq — > gg), as 
in Sec. 4.1. The resulting anomalous dimension matrices are given by Eqs. (^) and (|4^). We 
must diagonalize these matrices, finding eigenvalues Aq, and the transformation matrices Rj^, 
Eq. (p2|), which are used to take the hard and soft functions to the diagonal basis space, as in 
Eq. (0). Finally, the exponents Eaf3 that appear in the cross sections are found from Eq. (|2^). 

As an example, we consider predictions for high-py jets at zero hadronic cm. rapidity {1] = 0) 
in pp scattering at ^/s = 1.8 TeV. We work in valence quark approximation, and therefore need 
only the exponents for quark- ant iquark initial states. For Q, we take the patch discussed in Sec. 
4.1, whose exponents E^is were plotted in Figs. ^ and ^. The resulting ratio p is shown in Fig. 
D for pt = 50, 300 and 500 GeV as a function of x = Qn/pr, the ratio of the transverse energy 
flow to Pt, for < X < 0.4. (We set p to zero below Aqcd-) For larger Qq, we expect recoil. 
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Figure 5: The distributions p{x = Qn/pr) for the region Q defined in Sec. 4.1, with the observed 
jet at ?7 = with fip = pr = 50, 300 and 500 GeV. This calculation was carried out in valence 
quark approximation for pp scattering, at y/s = 1.8 TeV. 



which we have not taken into account, to be important. In this calculation, we have chosen the 
renormalization and factorization scales equal to pt, and have used CTEQ5L parton distribution 
functions . The dependence of the ratio p on the choice of the factorization scale is very weak. 

Recall that p{x) measures the fraction of events for which the radiated transverse energy is 
less than or equal to xpT- For any choice of fl the real parts of the exponents E, Eq. (^3]), are 
positive. As a result, the ratio p vanishes for Qq = Aqcd, and rises rapidly with x. For example, 
when Pt = 500 GeV, roughly 80% of events have Qq < 50 GeV for this choice of fl. In addition, 
for the choice of Q given above, lie Ea/B is almost always less than unity, and the distribution 
dp/dQn, found from Eqs. (|27|) and (|26|) , actually diverges at Qn = Aqcd- A similar behavior 
was found for "gap" dijet events, with low radiation in a large region that covers all < < 27r 
between high-p^ jets that are widely-separated in rapidity |jl5[. In such gap events there is an 
anomalously small eigenvalue, corresponding to a color exchange that approaches pure singlet 
in the t-channel for rj —>■ oo. In contrast, for the region Q considered here, which is outside the 
scattering plane and of limited extent in rapidity, all the eigenvalues are relatively small. 

The lowest-order expansion of the differential cross section dcr^^^ / dfjdQfi, Eq. corresponds 
to the eikonal approximation for the NLO correction to the cross section [|l^]. This cross section 
is proportional to the real parts of the exponents Eai3, Eq. (P^D, and is hence proportional to the 
size of the region Q in rj — (p space. The vanishing of p at low energy flow is an expression of 
the tendency for the scattered quarks to radiate, essentially according to the "antenna pattern" 
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Figure 6: The 77-dependence of p{x) for a fixed region Q defined as in Sec. 4.1., for radiated 
energies Qn = 1 GeV and 10 GeV, respectively, and pr = 300 GeV. 



TT] , |20| , Pq] of the corresponding classical field. The divergence of the perturbative one-loop 
running coupling is responsible for the vanishing of p{x) for x = Aqcd/pt rather than at zero, as 
would be the case for a fixed coupling. The rapid rise in p{x) indicates that most events generate 
relatively small amounts of perturbative radiation into the patch region. The differences in the 
curves are due to the running coupling, and become less important as pt increases. Finally, in 
Fig. H, we show the r^-dependence of p for two values of Qn for the same patch, with px = 300 
GeV. The rapidity-dependence is mild for this range in t] and choice of Q. 



5. Discussion 

The few examples we have given in the previous section illustrate the kind of new results on 
energy flow that can be derived using the methods outlined in this paper. In principle, given 
any fixed set of jets, or heavy quarks [^, of specified momenta and directions, we may calculate 
the distribution of perturbative energy flow into any fixed region Q of rapidity and azimuthal 
angle in the eikonal approximation. After resummation in logarithms of Qq/pt, the distribution 
is conveniently described in terms of the ratio piQu/pr), which specifies the percentage of dijet 
events for which a transverse energy of no more than Qq flows into region Q. 

The corresponding distribution, dp/drjdpTdQn may also be interpreted as a correlation of 
transverse energy flows, between the jet cone and region Q, generalizing energy-energy corre- 



lations to "moving cone" multijet algorithms, introduced by Giele and Glover |Q. These 



23 



are examples of Tkachov's "jet discriminators" |]3T[. In the formalism of pO|, the discussion in 



this paper relates to a three-cone configuration, in which two cones sample high-p^- jets, while 
one "cone", covering Q, encounters relatively soft radiation. We also anticipate that appropriate 
moments with respect to transverse energy in Q may be analyzed in the spirit of to gain 
insight into the influence of nonperturbative corrections for this class of correlations. 

In this paper, we have investigated only the relatively simple case of valence approximation. 
Beyond this, we hope that the kind of energy flow analysis described here will provide an addi- 



tional tool to study more complex jet events and production mechanisms at short distances |12 



Our predictions are at present limited by our neglect of "non-global" logarithms, which may be 
studied numerically [|171. We anticipate that by specifying restrictions on radiation into the in- 



termediate region Q more completely, it will be possible to regain a closed form for contributions 
of radiation into fl to certain inclusive cross sections. These and related issues are under study. 
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